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Abstract 

Consider an infinite tree with random degrees, i.i.d. over the sites, with a pre- 
scribed probabiUty distribution with generating function G{s). We consider the fol- 
lowing variation of Renyi's parking problem, alternatively called blocking RSA: at 
every vertex of the tree a particle (or "car") arrives with rate one. The particle sticks 
to the vertex whenever the vertex and all of its nearest neighbors are not occupied 
yet. 

We provide an explicit expression for the so-called parking constant in terms of 
the generating function. That is, the occupation probability, averaged over dynamics 
and the probability distribution of the random trees converges in the large-time limit 
to (l-a2)/2 with/„^^ = l. 
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1 Introduction 

In the classical car parking problem considered by Renyi ^llj, one dimensional cars with 
unit length appear one by one with their midpoints uniformly distributed over an interval. 
A car is parked unless it intersects with one or more previously placed cars. The process 
stops, when there is no further possibility of placing a car. In the discrete version of 
Renyi's car parking problem, cars of length 2 try to park at their midpoints randomly on 
the integers. Now a car can be parked if the distance of its midpoint to all other midpoints 
of already parked cars is two or more. It is well known that the probability that a given 
site is occupied by the midpoint of a car converges to (1 — e~^)/2 in the fully parked 
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state, that is when time tends to infinity [21 18]. In higher dimensions a model of this type 
will in general not be solvable and more complicated behavior is expected to occur. In 
physical language, this process is called (blocking) RSA (random sequential adsorption), 
the motivation being that particles are deposited onto a surface [21 [ini El [6]. For a nice 
review article over various packing problems (of which the parking problem can also be 
seen as a special example), see [1]. 

In this note we consider the problem of parking on a random tree, where particles 
at the vertices appear with an exponential waiting time, and parking is prohibited when 
there is a car-midpoint already at the site or at a neighboring site on the tree. We stress 
that the tree itself is unchanged under the dynamics. We solve the model and calculate 
the occupation probability as a function of time for random trees with i.i.d. degrees, 
averaged over the distribution of trees. In particular we recover the regular tree result 
[lOj . We believe that the random tree model has appeal for two reasons: The first is that 
it provides a step to the analysis of the process on random networks, motivated e.g. by 
the study of communication networks. Indeed, many of the random graph distributions 
proposed for its description which have found recent interest allow for a local random 
tree approximation [31 HI [7]. The second reason why it is nice is just the mathematical 
simplicity of the result which should not go unnoted. 

The proof is based on an analysis of the ordinary differential equation of the occupation 
probability at a fixed site of a fixed realization of the random tree. This is like [101 [8], 
but we give a self-contained exposition. Conditioning on a non-arrival of a car at this 
site the r.h.s. of this differential equation allows for a factorization in terms of occupation 
probabilities on rooted trees. Again a differential equation for these quantities can be 
derived. While it cannot be expected to be solvable for a particular realization, the 
averaging over the tree distribution still allows for closed-form expressions, as we shall see. 

1.1 Random Trees with Independent Identically Distributed Degrees 

The precise definition of the model is as follows. We consider a random tree with vertices 
i and degree at the site i given by Di. We choose Di to be independent random variables 
with the same distribution Q given by 

Q(A = k) = ak (1.1) 

on the integers starting from 2. The latter requirement ensures that we have no open ends 
with probability one. We denote the generating function of the distribution by 

oo 

G{s) = Y^akS^ (1.2) 

k=2 

We will denote the expected value with respect to this probability distribution by the 
same symbol Q. 
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1.2 The Dynamics 

Fix a tree with vertex set V. For any such tree we wih define a Markov jump process 
on the occupations numbers n = {ni)i^v € = {0, 1}^. Here denotes the occupation 
number of vertex i, meaning that 



if vertex i is vacant 

1 if vertex i is occupied 



The dynamics of the process is defined in terms of the generator which is given by the 
r.h.s. of the differential equation 



Jt 



rk[n) 



(1.3) 



with 



and 



rk{n) 



n,- 



rii \l k ^ i 
1 {ik = i 

- n <' 

i:d{i,k)<l 
= 1 - Ui, 



(1.4) 



and where d denotes the distance on the tree. 

This generator defines a Markov jump process on the infinite graph by standard theory 
[9], such that (|1.3p holds for any local function / : ^ M. Here E" denotes the expected 
value with respect to the process, started at the initial configuration n. 

This generator describes the parking of cars at all possible sites k at rate 1 given a 
configuration which is given by n. This parking is possible whenever the site and its 
nearest neighbors are vacant, that is rfc(n) = 1. Since at most one car can be parked at 
a given site, it suffices to consider the time of the first arrival of a car at this site, for all 
sites. 



1.3 Results 

We provide an explicit integral formula for the density of occupied sites at time t, averaged 
with respect to dynamics E and tree distribution Q. 



Theorem 1 



QE(no(t)) 



1 - a^(l 



where the function a{u) is defined by the equation 

xdx 



u 



(1.5) 



1.6) 



4u) G'(x) 

In particular, the occupation probability, averaged over dynamics and the probability distri- 
bution of the random trees converges in the large-time limit to (1 — a^)/2 with = 1. 
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E(no(t)) 



Remark: Let us specialize to the deterministic case of a regular tree with D > 2 
nearest neighbors. We obtain 

' i(i_e-2(i-e-*)) ifD = 2 

i(l - (1 + (Z) - 2)(1 - e-*))"^) if D > 2 
and recover in this way the known results on the integers and general regular trees [21 [HI [10] . 

2 Proof 

Consider a fixed realization of the random tree and look at the time-evolution given by 
(jl.3p . We use a short notation for the probability that a given set of sites A is empty at 
time t, 

a(A):=E(J]<(t)) (2.1) 

Similarly we write for the conditional probability that a given set of sites A is empty, 
conditional on the event that on another set of sites B no particles have arrived yet: 

Ct{A\B) := E([]<(n)|T,- > t for j G B) ^2.2) 

Of course, if no particle has arrived at a site it is empty, but the converse is not true. 
Here Tj denotes the arrival time of the particle at site j. The Tj's are i.i.d. exponentially 
distributed with expected value equal to one. We also write in short Ct{ii, ■ ■ ■ ,ik) = 
Ctiih, . . .,ik}) for A = {h, ...,ik} etc. 

As the underlying tree is random, also these correlation functions are random variables 
with respect to the distribution Q. For any realization we can write the dynamics ()1.3p as 

-j^CtiA) = Ct(A) (2.3) 

where A = {i £ G\d{i,A) < 1}. Similarly we have for the dynamics of the conditional 
correlation function the differential equation 

-±C.(A\B) = C.(1|B) p^^j 
= CtiAnB''\B) 

The first equation follows from restricting the generator to the set B'^ (since in B no 
particles have arrived yet). The second equation expresses the fact that when no particles 
have arrived in the set B, it implies that AnBis unoccupied. 

Fix an arbitrary vertex. Let us call this vertex 0. Starting from the dynamics for the 
correlation functions (12.31) we have 



■j^CtiO) = Cti{i,d{i,0)<l}) 

=e( n <(t)iToM)+E( n <(oiTo<t) ^^'^^ 

i,d(j,0)<l i,d{i,0)<l 



4 




-1 



Figure 1: Fragment of a tree where sites and 1 are empty and the occupancy of the other 
sites is unknown. 



We notice that the second term vanishes since an arrival of a particle at would imply 
that the sites {i,d{i,0) < 1} cannot all be empty. Writing the first term as a conditional 
probability we get 

-^Ct{0)=Ctai,d{i,0) < l}10)e-* (2.6) 
at 

Conditional on the event that no particle has arrived at the site 0, the dynamics for the 
occupation numbers on the branches of the tree that are emerging from is independent. 
Consequently the correlations factorize into a product and we have 

Ct{{i, d{i,0) <l}\0)= H Qm (2.7) 

i,d{i,0)=l 

where the Cj(i|0) is the non-occupation probability of a site i that is adjacent to the root 
on a rooted tree, assuming that no particle has arrived at the root. 

For a particular realization of the underlying tree, these functions will in general differ. 
Let us however now take the expected value over the tree distribution Q. We can decom- 
pose this average into an average over Dq (the number of neighbors of the central vertex 
0) and a conditional average over the remaining branches, rooted at 0. Conditioning the 
number of branches emerging from zero to be equal to k, the correlation functions are 
independent random variables with respect to Q and so we obtain 

Q(Ct{{i,d{i,0) = 1}\0)\Do = k^= {Q(Ct{l\0))y (2.8) 
where we denoted by the site 1 one of the nearest neighbors of 0, see Fig. [TJ 
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From here, ()2.6p and ()2.7p we have 



-^Q(a(0)) =Y,a,{QCtm) 

k=2 y^-^) 

= G{QCt{m)e-' 

Integrating this equation gives us the desired function Q(Cf(0)) once the "half-infinite 
function" QCt(l|0) is known. So, it remains to solve a problem on a rooted tree, condi- 
tioned that no particle has arrived at the root, averaged over Q. 

Indeed, it turns out to be possible to derive a closed-form differential equation for this 
object, as we show now. Note at first that it suffices to consider the event that no particle 
has arrived at 1 and write 

Ct{{i, d{i, 1) < 1}|0) = Ct{{i, d{i, 1) < 1}|0, l)e-* (2.10) 

Since the root and the site 1 itself are empty under the conditioning we have 

Ct{{i,d{i,l) < 1}|0,1) = Ct{{i + 0,d(i,l) = 1}|0,1) 

= n ^^-^^^ 

i7^0,d(i,l)=l 

Employing the dynamics (|2.4|) . using (|2.1U|) and averaging over Q we thus have 



-^QCi(l|0) =QCt({i,(i(i,l) < 1}|0) 



)( n ^t(^|0'l))e" 

i^O,d(i,l)=l 



t 



(2.12) 



Using conditional independence over the branches with respect to the tree distribution Q 
and noting that QCt(i|0, 1) = QCf(lIO) for the i's appearing in the above product to get 



J] Ct(i%\)) =Y,amCt{mf-^ (2.13) 

j7^0,d(i,l)=l k=2 

we finally obtain with (j2.12p and (j2.13p the closed-form differential equation 

-^Oafl|0)-^i^^™le- (2 14) 

JP-^'^^\^^ - Qa(l|0) ' ^ ^ ^ 



Its solution QCj(l|0) with initial condition QC(=o(l|0) = 1 is given by the integral 

1 - e"* (2.15) 



QCt{i|o) G{y) 
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Let us substitute this solution into the equation for the filhng probabihty (j2.9p . This 
then imphes the explicit formula for the Q-average of the fihing density of a random tree 
with generating function G of the form 

QE(no(t)) = f G{y{s))e-'ds (2.16) 
Jo 



where y{s) is defined by the equation 

xdx 

l-e^' (2.17) 



ly{s) G{x) 

Finally, with the substitution n = 1 — and recalling ()1.6p we see that 



G{y{s))e-'= / Giyiu))du = - / ydy 



^0 



(2.18) 

where the second step follows from ()2.17p . This finishes the proof. 
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